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Abstract 

We construct Hadamard matrices of orders 4 • 251 = 1004 and 4 • 631 = 2524, and skew- 
Hadamard matrices of orders 4 • 213 = 852 and 4 • 631 = 2524. As far as we know, 
such matrices have not been constructed previously. The constructions use the Goethals- 
Seidel array, suitable supplementary difference sets on a cyclic group and a new efficient 
matching algorithm based on hashing techniques. 

1 Introduction 

There are only 13 integers v < 500 for which no Hadamard matrices of order Av are known, see 
[U |2j [8]. In this paper we are able to remove the integer v = 251 from this list by constructing 
a Hadamard matrix of order 4 • 251 = 1004 by constructing cyclic supplementary difference sets 
(SDS) with parameters (251; 125, 120, 115, 115; 224). Hadamard matrices of order 4-631 = 2524 
are not known to exist, according to the table on page 278 of the second edition of the Handbook 
of Combinatorial Designs pQ , see also [U [6J . In this paper we construct Hadamard matrices of 
order 4 • 631 by constructing cyclic SDS with parameters (631; 315, 330, 330, 330; 674). 

Skew-Hadamard matrices of orders 4 • 213 and 4-631 are not known to exist either, see [TUj 
and [3] for the list of known orders of skew-Hadamard matrices. In this paper we construct 
skew-Hadamard matrices of orders 4-213 and 4-631 by constructing cyclic SDS with parameters 
(213; 106, 106, 105, 92; 196) and (631; 315, 330, 330, 330; 674), respectively. 

All SDSs in this paper were constructed by using the well-known method of taking base blocks 
to be unions of orbits of an automorphism group of the underlying cyclic group and a new effi- 
cient matching algorithm. This new algorithm turned out to be a crucial part of the discovery 
of these new matrices. 

A Hadamard matrix of order 428 has been constructed in [9] and currently the smallest multiple 
of 4 for which a Hadamard matrix is not known is 668. We refer the reader to [7] and [TUJ for 
more information on Hadamard and skew-Hadamard matrices. 
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The updated list of integers v < 500 for which no Hadamard matrices of order Av are known 
consists now of 12 integers 

167, 179, 223, 283, 311, 347, 359, 419, 443, 479, 487, 491 

all of them primes congruent to 3 (mod 4). 

2 The matching algorithm 

One of the serious difficulties of the unions of orbits approach when used in conjunction with 
the Goethals-Seidel array is that one needs to locate four lines in four text files such that the 
element- wise sums of these four lines are all equal to the same (pre-defined) constant. The 
difficulty stems from the fact that if two files contain ten million lines each, then the file with 
all possible sum combinations of these two files will contain 10 14 lines and therefore the ensuing 
naive algorithm to solve this problem, i.e. to combine the files two by two in two pairs and 
look for a potential match, is utterly impractical. In order to circumvent this difficulty, one 
potential way to avoid this combinatorial explosion is to resort to hashing functions. 

Let us begin with a precise statement of the problem we were faced with. Given four text files 
with M columns and Aq, N 2) N 3 , N 4 lines each, one has to identify four lines, one line in each 
file, identified by their line numbers in each file: li, l 2 , / 3 , U, such that their M element-wise 
sums are all equal to a constant A. In the case of the SDS (631; 315, 330, 330, 330; 674) it turns 
out that M = 21, A = 674 and that the four files contained approximately one million lines each. 

Let Z™ denote the additive monoid of non-negative integer n-tuples. Given four subsets 
A 1 ,A 2 ,A 3 ,A 4 G Z™ and a target s G Z™, we need to find a, L G A i: % — 1,2,3,4, such that 
0.1 + 02 + 03 + 04 = s, or prove that they do not exist. 

A naive exhaustive search requires 0(Y[ \Ai\) time and constant space. A meet-in-the- middle 
approach reduces this to 0(|Al| \A 2 \ + \A 3 \ \A 4 \) time but takes 0(|A:l| \A 2 \) space, by first storing 
all the sums a x + a 2 , where Oi G A ± , a 2 G A 2 , in a hash set H, and then searching for a 3 G A 3 , 
a 4 G A 4 , such that s — a 3 — a 4 G H . 

For convenience, we took B\ = { |_|J — a 3 | a 3 G A 3 } and B 2 = { |~|] — a 4 | a 4 G A 4 } and reduced 
the original problem to the equivalent one of finding a x G A x , a 2 G A 2 , b x G B 1 , b 2 G B 2 , such 
that di + a 2 = bi + b 2 . In our case, the elements of B\ and B 2 turned out to be non- negative 
integer tuples. 

In order to reduce space requirements and speed up additions, we mapped the n-tuples to 
64-bit non-negative integers using a linear hash function h : Z" — >■ Z/2 64 Z. Then we listed the 
solutions of h(ai) + h(a 2 ) = h(bi) + h(b 2 ), and for each of them checked whether they yield a 
solution to the original problem (some of them may not because of collisions). 
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Furthermore, we parallelized the algorithm, reducing the overall search time and the space 
requirements on each worker by a factor of M, where M is the number of workers. First, we 
represented each of the four input sets (after hashing) as an array of lists, which at index i, 
< i < M, stores the list of all elements equal to i modulo M. Using this data structure, 
worker number i, < i < M, can easily enumerate and store in the hash set H all sums a± + <22 
equal to i modulo M, and then lookup in H all sums b\ + 62 that are also equal to i modulo 



In practice, M was chosen to be much greater than the available number of workers, to further 
reduce memory use. Each worker pre-loaded the four sets in memory (represented with arrays 
of lists, as described above) and sequentially processed several remainders modulo M, where 
the number of remainders to be processed in a single run was taken large enough, so that the 
search time would dominate the pre-loading time, but would not exceed the maximal allowed 
duration of a job on the cluster. Reduced memory requirements enabled us to simultaneously 
schedule multiple workers on a single multi-core machine, thus fully utilizing the capacity of 
the cluster. 

3 Results 

We now present our results in the form of SDSs, for v = 213, 251, 631. Non-equivalence of SDSs 
was established by an implementation of the method described in [5]. 

We define the notation for the orbits of the action of the subgroup that we use to construct 
the solutions below. The automorphism group of the additive cyclic group Z v will be identified 
with the group of invertible elements, Z*, of the ring Z v . Consider a fixed subgroup H of order 
h of Z*. Clearly h must divide \Z*\ = <fr(v). Denote by H ■ k the orbit of H in Z v through 
the point k, where • is multiplication mod v. We refer to the orbit H ■ = {0} as the trivial 
orbit. The orbit H ■ 1 is just the subgroup H itself. In general the size of an orbit may be any 
divisor of \H\ and if v is a prime then every nonzero orbit is just a coset of H in Z* and so will 
have size \H\. 
The notation 



will be used below to present all the solutions found, in fact each solution will be given only via 
the four index sets J, K, L, M. For suitable choices of the four index sets J, K, L, M, the four 
sets X, Y, Z, W defined in ([T]) form SDS(i>; x, y, z, w; A) with x = ^,- 6 j \H-j\,y = J2keK \H-k\, 
z = ^2 leL \H ■ l\, w = J2meM \H ' m \ an d X = x + y + z + w-v. 



M. 




(1) 
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All Hadamard matrices in this paper are constructed via the Goethals-Seidel array 



Pi 


P 2 R 


P3R 


P 4 R 


P 2 R 


Pi 


-P?R 


P3R 


P 3 R 


PJR 


Pi 


-P?R 


P 4 R 


-P?R 


P?R 


Pi 



where R denotes the v xv matrix with ones in the back-diagonal and zeros everywhere else. To 
obtain a Hadamard matrix via the Goethals-Seidel array one has to substitute the v x v matrices 
Pi,P 2 , P3, -P4 by the {±1} circulant matrices that arise from the four subsets X, Y, Z, W that 
make up the SDS(t> ; x, y, z, w; A). More precisely, let ax = (a , . . . , a„_i) be a sequence defined 
by a,i = — 1 if i G X, and = 1 if % X; and define the sequences ay,az, aw similarly Then 
denote by [X], [Y], [Z], [W] the {±1} circulant matrices whose first rows are ax, ay, az, aw 
respectively. If X, Y, Z, W form a SDS(i>; x, y, z, w; A), then 

[X][X] T + [Y][Y] T + [Z][Z] T + [W][W] T = 4vl v . 



3.1 ^ = 213 

Consider the subgroup H = {1, 37, 91, 103, 172, 187, 190} of order 7, of Z* 13 . Note that Z* 13 is 
of order 0(213) = 140. Here are the 32 nontrivial orbits of the action of H on Z213. 



H 


1 = 


[1,37,91,103,172,187, 190] 


H 


21 


= [21,33,93,138,156,204,207] 


H 


2 = 


[2,74,131,161,167, 182,206] 


H 


22 


= [22,67,85,133,136,163,175] 


H 


3 = 


[3,60,90,96,111,135,144] 


H 


23 


= [23,26,41,110,122,176,212] 


H 


4 = 


[4,49,109,121,148,151,199] 


H 


28 


= [28,115,124,130,184,205,208] 


H 


5 = 


[5,8,29,83,89,98,185] 


H 


30 


= [30,45,48,72,108,162,174] 


H 


6 = 


[6,9,57, 75,120,180,192] 


H 


34 


= [34,70,94,97,112,181,193] 


H 


7 = 


[7,31,46,52,82,139,211] 


H 


38 


= [38,50,77,80,128,146,191] 


H 


10 = 


[10,16,58,157,166,178,196] 


H 


39 


= [39,51,105,141,165,168,183] 


H 


11 = 


[11,68,140,149,173,188,194] 


H 


42 


= [42,63,66,99,186,195,201] 


H 


12 = 


[12,18,27,114, 147, 150,171] 


H 


43 


= [43, 76, 79, 100, 154, 160, 169] 


H 


13 = 


[13,55,61,88,106,118,127] 


H 


44 


= [44,53,59,113,134,137,170] 


H 


14 = 


[14,62,65,92,104,164,209] 


H 


69 


= [69,78,102,117,123,153,210] 


H 


15 = 


[15,24,36,54,81,87,129] 


H 


71 


= [71] 


H 


17 = 


[17,35,47,56,155,197,203] 


H 


84 


= [84, 126, 132, 159, 177, 189, 198] 


H 


19 = 


[19,25,40,64,73,145,202] 


H 


86 


= [86,95,107,125,152,158,200] 


H 


20 = 


[20,32,101,116,119,143,179] 


H 


142 


= [142] 



We give an SDS with parameters (213; 106, 106, 105, 92; 196), via the index sets J, K, L, M (with 
respective cardinalities 16,15,14,16) to be used in (JI|), which gives rise to skew-Hadamard 
matrices of order 4 • 213 = 852. The associated decomposition into the sum of 4 squares is: 

4 • 213 = 852 = l 2 + l 2 + 3 2 + 29 2 = 

= (213 - 2 ■ 106) 2 + (213 - 2 ■ 106) 2 + (213 - 2 • 105) 2 + (213 - 2 ■ 92) 2 . 
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J = {4,5,7,10,11,15,17,19,20,30,34,38,39,42,43,142} 

K = {2,7,12,14,17,22,28,34,39,42,43,44,69,84,86} 

L = {1,4,11,15,17,20,21,28,30,34,42,44,69,142} 

M = {2,4,10,11,12,15,21,22,23,28,30,34,44,69,71,86} 



3.2 v = 251 



Consider the subgroup H = {1, 20, 113, 149, 219} of order 5, of Z% 51 . Note that Z\ hX is of order 
0(251) = 250. Here are the 50 nontrivial orbits of the action of H on Z251. 



T T 

H 


1 


= 


T"i c\r\ -1-10 -1 a r\ ai a! 

[1, 20, 113, 149, 219] 


T T 

H 


42 


= 


[42, 87, 162, 228, 234] 


H 


32 


= 


[32, 102, 138, 231, 250] 


H 


18 


= 


[18, 26, 109, 172, 177] 


T T 

H 


2 


— 


[2, 40, 47, 187, 226J 


T T 

H 


74 


— 


[74, 79, 142, 225, 233] 


J T 

H 


25 




for n a c\(\ \ 011 c\ a c\\ 

[25, 64, 204, 211, 249] 


J J 

H 


19 




l"i c\ f~7 r\ 1 on 1 on 1 a f~\ 

[19, 70, 129, 139, 145] 


J T 

H 


3 




[3, 60, 88, 155, 196J 


J T 

H 


106 




[106, 112, 122, 181, 232] 


T T 

H 


55 


= 


r f f /"\ c 1 /u) 1 c\i r\ a c~>~\ 

[55, 96, 163, 191, 248J 


H 


21 


= 


[21, 81, 114, 117, 169] 


H 


4 


— 


[4,80,94, 123,201] 


H 


82 


= 


[82, 134, 137, 170,230] 


T T 

H 


50 


— 


r f r\ 1 c\c~> 1 f ^7 1^1 c\ A ^rl 

[50, 128, 157, 171, 247] 


T T 

H 


28 


— 


[c\ f 1 r\ 1 f c\ I p*" /"* 1 

[28, 58, 108, 152, 156] 


T T 

H 


5 


= 


r f /"> f "> r\ 1 -1 r\r\ c\ a ol 

[5, 63, 91, 100, 243] 


T T 

H 


C\ f 

95 


— 


\ r\ f /\ ( \ 1 /io ir\o nool 

[95, 99, 143, 193, 223] 


H 


8 


= 


[8, 151, 160, 188, 246] 


H 


30 


— 


[on a a f\r~i I c\'t onnl 

[30, 44, 98, 127, 203] 


T T 

H 


r* 






\ n rn 1 on 1 a 1 1 l ~7r i ~\ 

[6, 59, 120, 141, 176] 


T T 

H 


A O 

48 




[48, 124, 153, 207, 221] 


H 


75 




[75,110,131,192,245] 


H 


33 




[33,148,158,199,215] 


H 


7 




[7, 27, 38, 39, 140] 


H 


36 




[36,52,93,103,218] 


H 


111 




[111,212,213,224,244] 


H 


34 




[34,46,77,167,178] 


H 


9 




[9,13,86,180,214] 


H 


73 




[73,84,174,205,217] 


H 


37 




[37, 71,165,238,242] 


H 


35 




[35, 135, 190, 195, 198] 


H 


10 




[10,126,182,200,235] 


H 


53 




[53,56,61,116,216] 


H 


16 




[16,51,69,125,241] 


H 


41 




[41,67,85,115,194] 


H 


11 




[11,133,150,220,239] 


H 


57 




[57,136,166,184,210] 


H 


12 




[12,31,101,118,240] 


H 


43 




[43,90,107, 130,132] 


H 


14 




[14, 29, 54, 76, 78] 


H 


119 




[119,121,144, 161,208] 


H 


173 




[173,175,197,222,237] 


H 


45 




[45, 65, 66, 147, 179] 


H 


15 




[15,22,49,189,227] 


H 


72 




[72, 104, 185, 186, 206] 


H 


24 




[24, 62, 202, 229, 236] 


H 


68 




[68,83,92,105,154] 


H 


17 




[17,23,89,164,209] 


H 


97 




[97, 146, 159, 168, 183] 



We give two SDSs with parameters (251; 125, 120, 115, 115; 224), via the index sets J, K, L, M 
(with respective cardinalities 25,24,23,23) to be used in ([1]), which give rise to Hadamard 
matrices of order 4 • 251 = 1004. The associated decomposition into the sum of 4 squares is: 



4-251 = 1004 



1 2 + 11 2 + 21 2 + 21 2 = 
(251 - 2 • 125) 2 + (251 



2 • 120) 2 + (251 - 2 • 115) 2 + (251 - 2 • 115) 2 . 
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J = {2, 4, 5, 6, 7, 111, 9, 37, 10, 11, 12, 14, 173, 15, 24, 74, 19, 48, 33, 73, 53, 57, 43, 72, 68} 

K = {2, 55, 6, 75, 16, 11, 173, 24, 17, 42, 18, 19, 106, 21, 30, 48, 34, 73, 35, 53, 41, 57, 43, 97} 

L = {1,2,4,50,6,75,9,37,10,11,173,15,17,18,21,48,36,34,119,45,72,68,97} 

M = {32, 2, 25, 3, 5, 8, 7, 16, 173, 15, 74, 19, 106, 21, 95, 30, 73, 53, 41, 57, 45, 72, 68} 

J = {32, 3, 55, 50, 5, 6, 7, 9, 37, 11, 12, 14, 15, 24, 18, 19, 82, 30, 36, 34, 73, 35, 41, 57, 72} 

K = {1, 32, 25, 55, 5, 8, 75, 9, 37, 11, 14, 15, 24, 17, 42, 18, 19, 95, 48, 41, 57, 43, 45, 97} 

L ={1,2, 25, 3, 50, 8, 75, 9, 11, 12, 173, 24, 17, 28, 95, 30, 33, 36, 34, 73, 45, 72, 97} 

M = {32, 25, 4, 5, 6, 75, 7, 16, 11, 12, 24, 42, 106, 28, 95, 48, 73, 35, 53, 41, 57, 45, 72} 

3.3 v = 631 

Consider the subgroup ff = {1, 8, 43, 64, 79, 188, 228, 242, 279, 310, 339, 344, 512, 562, 587} of or- 
der 15, of Z£ 3l . Here are the 42 nontrivial orbits of the action of if on Z^i. 

ff-1 = {1,8,43,64,79,188,228,242,279,310,339,344,512,562,587} 

if -2 = {2,16,47,57,86,128,158,376,393,456,484,493,543,558,620} 

if -3 = {3,24,53,95,129,192,206,237,274,299,386,401,424,499,564} 

ff -4 = {4,32,94,114,121,155,172,256,281,316,337,355,455,485,609} 

if. 5 = {5,36,40,133,215,286,288,309,320,395,411,433,458,509,579} 

H ■ 6 = {6, 48, 106, 141, 171, 190, 217, 258, 367, 384, 412, 474, 497, 548, 598} 

H-7 = {7,54,56,60,148,277,301,323,334,429,432,448,480,515,553} 

ff-9 = {9,10,72,80,159,191,235,266,285,387,430,527,572,576,618} 

ff ■ 11 = {11, 73, 88, 138, 147, 175, 238, 255, 473, 503, 545, 574, 584, 615, 629} 

ff ■ 12 = {12, 96, 103, 137, 193, 212, 282, 317, 342, 363, 380, 434, 465, 516, 565} 

ff ■ 13 = {13, 55, 59, 104, 201, 244, 346, 365, 396, 440, 472, 551, 559, 621, 622} 

ff ■ 14 = {14, 15, 37, 108, 112, 120, 227, 233, 265, 296, 329, 399, 475, 554, 602} 

if -17 = {17,41,81,84,89,90,100,136,169,222,326,328,457,501,514} 

if • 18 = {18, 20, 143, 144, 160, 229, 318, 382, 423, 470, 513, 521, 532, 570, 605} 

if • 19 = {19, 131, 152, 181, 186, 211, 226, 239, 253, 263, 417, 426, 546, 582, 585} 

ff • 21 = {21, 25, 34, 82, 162, 168, 178, 180, 200, 272, 283, 338, 371, 397, 444} 

if • 22 = {22, 146, 176, 276, 294, 315, 350, 375, 459, 476, 510, 517, 537, 599, 627} 

ff • 23 = {23, 107, 184, 189, 196, 210, 225, 250, 306, 340, 358, 418, 518, 538, 555} 

ff -26 = {26,61,99,110,118,161,208,249,313,402,471,487,488,611,613} 

ff • 27 = {27, 28, 30, 74, 167, 216, 224, 240, 319, 454, 466, 477, 530, 573, 592} 

ff • 29 = {29, 77, 156, 232, 302, 335, 366, 398, 404, 511, 519, 523, 594, 616, 617} 

ff • 31 = {31, 71, 91, 97, 127, 145, 149, 248, 385, 413, 446, 529, 556, 561, 568} 

ff • 33 = {33, 83, 134, 157, 219, 247, 264, 373, 414, 441, 460, 490, 525, 583, 625} 

H ■ 35 = {35, 51, 109, 123, 241, 243, 252, 267, 270, 280, 300, 347, 353, 408, 507} 

ff • 38 = {38, 203, 221, 262, 304, 362, 372, 422, 452, 461, 478, 506, 526, 533, 539} 

ff • 39 = {39, 58, 101, 154, 165, 177, 312, 391, 407, 415, 464, 557, 601, 603, 604} 

ff • 42 = {42, 45, 50, 68, 111, 163, 164, 257, 324, 336, 356, 360, 400, 544, 566} 

H ■ 44 = {44, 69, 119, 287, 292, 321, 352, 389, 403, 443, 552, 567, 588, 623, 630} 

ff • 46 = {46, 49, 85, 205, 214, 368, 378, 392, 405, 420, 445, 450, 479, 500, 612} 
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H ■ 52 = {52, 122, 173, 198, 220, 236, 311, 322, 343, 345, 416, 498, 591, 595, 626} 

H ■ 62 = {62, 139, 142, 182, 194, 195, 254, 261, 290, 298, 427, 481, 491, 496, 505} 

H ■ 63 = {63, 70, 75, 102, 185, 218, 246, 383, 482, 486, 504, 534, 540, 560, 600} 

H ■ 65 = {65, 87, 231, 271, 275, 295, 307, 374, 467, 468, 520, 563, 581, 586, 589} 

H ■ 66 = {66, 115, 166, 197, 251, 268, 289, 314, 349, 419, 438, 494, 528, 535, 619} 

H ■ 67 = {67, 132, 207, 230, 245, 332, 357, 394, 425, 439, 502, 536, 578, 607, 628} 

H ■ 76 = {76, 93, 113, 213, 273, 291, 325, 381, 406, 421, 435, 442, 447, 524, 608} 

H ■ 78 = {78, 116, 151, 183, 199, 202, 297, 308, 330, 354, 483, 571, 575, 577, 624} 

H ■ 92 = {92, 98, 105, 125, 153, 170, 179, 209, 259, 269, 327, 369, 410, 428, 593} 

H • 117 = {117, 130, 174, 303, 305, 409, 462, 495, 531, 541, 542, 547, 550, 590, 614} 

H ■ 124 = {124, 223, 278, 284, 331, 351, 361, 364, 379, 388, 390, 508, 522, 580, 596} 

H ■ 126 = {126, 135, 140, 150, 204, 333, 341, 370, 377, 436, 437, 449, 489, 492, 569} 

H ■ 187 = {187, 234, 260, 293, 348, 359, 431, 451, 453, 463, 469, 549, 597, 606, 610} 

We give four nonequivalent SDSs with parameters (631; 315, 330, 330, 330; 674), via their index 
sets J,K,L,M (with respective cardinalities 21,22,22,22) to be used in ([[]), which give rise 
to Hadamard matrices of order 4 • 631 = 2524. In addition, the first two SDSs give rise to 
skew-Hadamard matrices of order 4 • 631 = 2524. The associated decomposition into the sum 
of 4 squares is: 

4 ■ 631 = 2524 = l 2 + 29 2 + 29 2 + 29 2 = 

= (631 - 2 ■ 315) 2 + (631 - 2 ■ 330) 2 + (631 - 2 ■ 330) 2 + (631 - 2 • 330) 2 . 



J = {1,2,3,4,6,7,12,13,14,17,19,21,26,27,31,38,42,52,62,76,124} 

K = {1, 2, 3, 4, 5, 6, 7, 9, 12, 17, 23, 26, 27, 31, 33, 38, 42, 44, 52, 76, 78, 126} 

L = {1, 2, 3, 4, 6, 7, 9, 11, 12, 13, 14, 17, 18, 19, 21, 29, 35, 46, 52, 62, 66, 76} 

M = {1, 2, 3, 4, 5, 6, 12, 13, 14, 18, 19, 21, 22, 26, 27, 38, 39, 42, 63, 67, 92, 124} 

J = {11, 13, 19, 22, 26, 29, 31, 33, 38, 39, 44, 52, 62, 65, 66, 67, 76, 78, 117, 124, 187} 

K = {1,3,4,5,9,11,14,17,18,22,23,26,29,33,38,39,42,46,62,65,67,117} 

L = {1, 2, 5, 6, 9, 17, 18, 21, 22, 27, 33, 39, 44, 46, 52, 66, 76, 78, 92, 117, 124, 187} 

M = {2, 4, 5, 7, 9, 11, 12, 13, 18, 19, 21, 23, 29, 31, 42, 44, 65, 66, 67, 78, 92, 187} 

J = {1, 3, 4, 7, 12, 13, 17, 18, 19, 27, 29, 31, 33, 35, 42, 46, 62, 67, 92, 124, 187} 

K = {5, 6, 7, 9, 12, 17, 18, 21, 22, 26, 29, 33, 35, 38, 42, 44, 62, 63, 67, 76, 124, 126} 

L = {4, 6, 7, 9, 11, 13, 14, 17, 18, 21, 22, 23, 27, 29, 35, 39, 46, 67, 78, 124, 126, 187} 

M = {2, 5, 6, 7, 11, 13, 14, 18, 19, 21, 27, 29, 33, 35, 39, 44, 52, 62, 63, 65, 117, 187} 

J = {1, 2, 4, 5, 6, 9, 11, 12, 13, 17, 22, 27, 29, 42, 44, 65, 78, 92, 117, 126, 187} 

K = {1, 7, 9, 13, 17, 19, 22, 23, 27, 29, 31, 33, 35, 38, 44, 46, 65, 66, 76, 78, 126, 187} 

L = {2, 3, 4, 5, 13, 17, 21, 22, 29, 35, 38, 39, 52, 62, 63, 65, 67, 76, 92, 117, 124, 126} 

M = {2, 3, 4, 5, 7, 9, 11, 12, 17, 18, 21, 22, 29, 33, 42, 44, 46, 52, 65, 66, 92, 187} 
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